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Abstract. Let ft C M", n > 3, be a smooth bounded domain and consider a 
coupled system in Q consisting of a conductivity equation V • "f{x)'Vu{t, x) = 
and an anisotropic heat equation Hi^^{x)dtip{t^x) = V • {A{x)Vi){t,x)) + 
(7Vu(t, a;)) • Vu{t,x), t > 0. It is shown that the coefficients 7, k and 
A = (ajk) are uniquely determined from the knowledge of the boundary map 
u\dn ^ I' ■ AViplgn, where is the unit outer normal to dQ. 

The coupled system models the following physical phenomenon. Given a 
fixed voltage distribution, maintained on the boundary dQ, an electric cur- 
rent distribution appears inside fl. The current in turn acts as a source of 
heat inside fl, and the heat flows out of the body through the boundary. The 
boundary measurements above then correspond to the map taking a voltage 
distribution on the boundary to the resulting heat flow through the bound- 
ary. The presented mathematical results suggest a new hybrid diffuse imaging 
modality combining electrical prospecting and heat transfer-based probing. 

Keywords: electrical impedance tomography, heat transfer, inverse problem, 
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1. Introduction 

Let us model a physical body by a bounded set Q C M", n > 3, with smooth 
boundary dQ, and the following spatially varying quantities: heat capacity c{x), 
density p{x), electric conductivity 7(x), and (possibly anisotropic) thermal con- 
ductivity A{x) = (ajfc(x)), each defined for x G fi. 

Consider applying a spatially and temporally variable electrical voltage distri- 
bution f{t,x) at the boundary dQ starting at time t = 0. Then, if there are no 
sinks or sources of current inside Q, the electric potential u{t,x) inside the body 
satisfies the conductivity equation 

V ■ -fVu(t, x) = for X G f2 and t > 0, 

~ d.l) 
u{t,.)\an = f{t,-). ^ ^ 

Equation (II .ip is often used as a mathematical model for electrical impedance 
tomography (EIT), where one measures the current through the boundary caused 
by a family of static voltage distributions f{t,x) = and recovers 7(x) from 
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such voltage-to-current map : (p ^ v ■ Vu\qq_. Here v is the unit outer 
normal to dVL. We refer to [22] for an extensive survey of the mathematical 
developments in EIT. See also [31 HJ [321 [361 EZ] for results in the two-dimensional 
case, and [I3l [271 |28l [291 EH [SI [38] for results in higher dimensional cases. For 
counterexamples to uniqueness of time-harmonic inverse problems involving very 
anisotropic and degenerate material parameters, leading to the phenomenon of 
invisibility, see [TOl [TTl [T2] . 

Our aim here is somewhat different as we wish to couple heat conduction to 
the problem. Let us denote the electrical power density inside Vt by F: 

F{t,x) = {-fVu{t,x)) ■Vu{t,x). (1.2) 

Now F acts as a source of heat inside Q. Assuming that the body is at a constant 
(zero) temperature at the time t = when the voltage is first applied, and the 
surface of the body is kept at that temperature at all times, the temperature 
distribution ^(t, x) inside Q satisfies the following heat equation: 

K-^{x)dtip{t, x) = V ■ (A(a;)V^(t, x)) + F{t, x) for x G and t > 0, ^ 

V'li+xan = 0, ^\t=o = 0, 

where k{x) = c{x)~^ p{x)~^ . The model (11.11) . (II. 2p . (II. 3p is based on the physical 
assumption that the heat transfer is so slow that the quasistatic (DC) model for 
the electric potential (11.11) is realistic. 

Associated to the coupled system (11.11) . (11.21) . (II. 3p . we introduce the voltage- 
to-heat flow map ^j,K,A defined by 

^^,^,A : f ^ 1^ ■ AV^l^^^g^. (1.4) 

The idea is to measure the heat flow through the boundary caused by the heat 
from the electric current resulting from the applied voltage distribution. 

Our main result is Theorem 12.11 below, stating that under certain smooth- 
ness assumptions, the coefficients 7, k, and A are uniquely determined from the 
knowledge of the voltage-to-heat flow map S^^^.a- 

The method of proof of Theorem 12 . 1 1 also outlines a constructive reconstruction 
procedure for recovering conductivity 7 from ^.A- Namely, it turns out that 
applying a temporally static voltage distribution f(t,x) = (t>{x) and studying 
^7,k,a/ at thermal equilibrium [t — )■ 00) yields the knowledge of the Dirichlet-to- 
Neumann map A^0 related to the EIT problem. Then one can recover 7 using 
Nachman's reconstruction result [31] . 

Notice that various hybrid imaging methods have been proposed and analyzed 
recently. Examples include thermoacoustic and photoacoustic imaging [H [SI [25] , 
combination of electrical and magnetic probing [221 [S3], electrical and acoustic 
imaging ^15] and magnetic and acoustic imaging [301 [2] . Theorem 12.11 suggests 
a new hybrid imaging method, utilizing two diffuse modes of propagation: elec- 
trical prospecting and heat transfer-based probing. We emphasize that the pro- 
posed method recovers complementary information about three different physical 
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Electrode for applying voltage 



Heat flow sensor 



Figure 1. Schematic illustration of the practical measurement 
setup motivating the proposed hybrid imaging method. The ideal 
voltage-to-heat flow map niay be approximated in practice 

by maintaining fixed voltages at the electrodes while measuring 
heat flow using the interlaced sensors. 

properties. We also note that in many applications where one wants to recon- 
struct the heat transfer parameters k{x) and A{x), the use of electric bound- 
ary sources may be easier than controlling the temperature or the heat flux at 
the boundary. Concerning inverse problems for the heat equation, we refer to 



We remark that in practice one might use a measurement setup shown in Figure 
[U However, analysis of such discrete measurements is outside the scope of this 
paper, and in the mathematical results below we work with the continuum models 



This paper is organized as follows. In Section |5] we state our assumptions and 
results in a mathematically precise form. In Section Owe give an auxiliary density 
result for the conductivity equation. Section H] is devoted to the reconstruction of 
the conductivity 7. The proof of Theorem 12. II is completed in Section [5l where we 
show the identiflability of the heat parameters k and A. Finally, Appendix A is 
devoted to the recovery of the boundary values of the matrix A from interior-to- 
boundary measurements, associated to a suitable elliptic boundary value problem. 
This result may be of an independent interest. 



Let C M", n > 3, be a bounded domain with C°° boundary. Let 7 G C°°{^1) 
be a strictly positive function on Q. Then given f{t,x) G C^(R^, H'^{dQ)), s > 
1/2, on the boundary at time t > 0, there exists a unique u G C^(]R+, 
which solves the boundary value problem 



I61CI[I91[201EIIE21[231S0]. 



dLl]), (OD, (O, and (m. 



2. Statement of results 



V ■ 7VM(t,x) = 
u{t,-)\9n = f{t, ■) 



for X G f2 and t > 0, 



(2.1) 
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see [16]. We have 

F{t, x) = {-fVu{t, x)) ■ Vu{t, x) e C\W+, L\n)), (2.2) 

provided that s is taken large enough, say s > (n + l)/2. In what follows we 
shall always choose the Sobolev index s in this way. Consider the anisotropic 
heat equation 

K'^{x)dt^/j{t, x) = V ■ (A(x)V^(t, x)) + F{t, x) iov xen and t > 0, 
^\w+xdn = 0, ^\t=o = 0. 

Here A{x) = {ajk{x)) is a real symmetric n x n matrix with aj±{x) G C°°{Q), 
and there exists Co > such that 

n 

J2 ('A^Mk > Coler, for all {x, G n X M^ (2.4) 

j,k=i 

We shall assume that < k G C°°{Q). The operator 

Pv = -k{x)V ■ {A{x)Vv) 

is formally self-adjoint in = L'^{Q, n^^dx) and we have 

{Pv,v)li>Co\\v\\1„ veC^iQ), Co>0. 

We also let P denote the Friedrichs extension of the operator P on C^{Q), so 
that the domain of the positive self-adjoint operator P is {Hq fl H'^){Q). 
The solution of (12.31) is given by the Duhamel formula 

ilj{t,x)= [ {e-^'~'^^KF){s,x)ds e C\R+,L\Q))nC{W+,{H^nH^){Q)), (2.5) 
Jo 

see [T6] . 

Associated to the coupled system (12. ip . (12.21) . and (12. 3p . we consider the 
voltage-to-heat flow map, 

The main result of the paper is as follows. 

Theorem 2.1. Assume that < 7^ G C°°(fi), < kj e C°°(fi), and Aj are real 
symmetric nx n matrices with C°°{Q) entries, satisfying (I2.4p . forj = 1,2. // 
= S^2,«2,A2; ^^en 7i = 72, Ki = K2 and Ax = A^. 

It turns out that in the course of the proof of Theorem 12. H we establish a 
result for the anisotropic heat equation, which may be of independent interest. 
In order to state the result, consider the inhomogeneous initial boundary value 
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problem fl2.3p for the anisotropic heat equation with an arbitrary source F G 
Ci(l+,L2(fi)). Define the map, 

F ^u- AV^)^^ 

where ip^ is the solution of fl2.3p . 



Theorem 2.2. Let C M", n > 2, he a hounded domain with C°° houndary. 
Assume that < kj E C°°(Q), and Aj are real symmetric n x n matrices with 
C°°{Q) entries, satisfying fl2.4l) . for j = 1,2. If 'E.^^^Ai = '^k2,A2j then Hi = ^2 
and Ai = A^. 

3. An auxiliary density result 

Let 17 C M", n > 3, be a bounded domain with C°° boundary and let 7 G 
C°°(i7) be a strictly positive function on We shall need the following density 
result, which is a quite straightforward consequence of 



Proposition 3.1. The set 

span{7Vwi ■ Vwa : Wj G C°^(n), V ■ 7Vwj = 0, j = 1, 2} 
is dense in L'^{Q). 
Proof Let / G L'^{Q) be such that 

f-fVwi-Vw2dx = 0, (3.1) 



n 

for any solution 1^1,^2 G C°°{Q) of the conductivity equation 

V-7Vu7 = 0. (3.2) 

We have 

7VW1 ■ Vw2 = -V ■ {'jV{wiW2)). (3.3) 

It follows from [38], see also [5j, that for any p G satisfying p ■ p = and 
IpI > 1 large enough, the conductivity equation (13. 2p has a solution 

w,{x)=e'P--^-'/\l + r^), (3.4) 

where Vp G C°°(fi) satisfies 

C 

\\rp\\H^in) < -TT, m > 0. (3.5) 

Here the constant Cm depends on Q, n, and a finite number of derivatives of 7. 

Given ^ G M" and R > 0, according to [38], there exist pi, P2 G C" such that 
Pj ■ Pi = 0, pi + P2 = ^ and \pj\ > R, j = 1, 2. 
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For the solutions Wp-^ and Wp^ of the form (13 ■4p . we get 

V ■ {lV{Wp,Wp,)) = V ■ (7V(7-'e'«-^(l + r,, + r,, + r,,r,J)) 

= (V ■ (7V7-^) + 7(V7-') • li - |ene^^-^(l + r,, + r,, + r,,r,J 
+ (7V7"^ + 2tOe'^-^ ■ {Vrp, + Vr,, + V(r,,r,J) 
+ e^«-^(Ar,,+Ar,,+A(r,,r,J). 

In view of (I3.3p . we may substitute the latter expression into (13. ip and let i? — )■ oo. 
Using (13. 5p . we obtain that 

/ (V ■ (7V7"') + 7(V7-^) ■ ii - \i\^)e'^-''fdx = for all ^ G M". (3.6) 
Jn 

Here we shall view 7 as a strictly positive C°° function on M", which is equal to 
a positive constant near infinity. The identity (13. 6p is equivalent to 

J-.^g(V ■ (7V7-')xq/ - V ■ (7V(7-')xn/) + A(xn/)) = 0, 

where J^x^^ denotes the Fourier transformation and xn is the characteristic func- 
tion of Q. It follows that xnf is a solution of a second order elliptic equation on 
M" with smooth coefficients. Since it is compactly supported, by unique contin- 
uation we conclude that / = in f2. This completes the proof. □ 

4. Recovering the conductivity 7 from the voltace-to-heat flow 

MAP 

The purpose of this section is to make the first step in the proof of Theorem 
12. H by establishing the following result. Recall that here n > 3. 

Proposition 4.1. The voltage-to-heat flow map ^■y,K,A determines the conduc- 
tivity 7 uniquely. 

When proving Proposition 14. ![ we let a G C°°(]R+; [0, 1]) be such that 

a|t<i/2 = 0, a\t>i = 1. 

Then set f{t, x) = a{t)h{x) with h G H'^{dQ), s large enough. Using the Duhamel 
formula (12. 5p . we shall study the behavior of ipit, x) as t +00. The solution u 
of (12. ip satisfies 

u{t, x) = a(t)wo{x), 

where Wq solves 

V ■ -fVwo(x) = in n, 
wo\an = h[x). 

Thus, 

F{t,x) = -fa^{t)Vwo ■ Vwo, 
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and (12 .Sp gives 

^(t,x)= [ a\s)e-'-'~'^^-fKVwo-Vwods. (4.2) 
Jo 

Lemma 4.2. We have 

\\ip{t, ■) - p-\-fnVwo ■ Vwo)IId(p) ^0, ast ^ +oo. 
HereV{P) = (if^nifo)(i7) is equipped with the graph norm \\ip\\v{P) = ||^||L2(n) + 

Proof. We shall first check that 

||^(t, ■) - P-\^kVwo ■ Vwo) |U2(f,) ^0, as t ^ +oo. (4.3) 
It follows from (14. 2 p that for t > 1, we have 

ilj{t,x) = / a^(s)e"^*~''^^7KVwo ■ Vwqc^s + / e'^^^nVwo ■ Vwods. 
Jo Jo 

Using that 

P + OO 

e-'^dt = p-\ 



in the sense of bounded operators on L^(fi), see [16], we get with L^-convergence, 
as t — +00, 

/ e '^^'^kVwq ■ Vwods — )■ P ^{^kVwq ■ Vi^o)- 
Jo 



On the other hand 

< I a'{s)\\e-^'-'^^-iKVwo-Vw^\\L2^^)ds 
Jo 

< / a^{s)\\e'^^^''^^\\c{L^n),L^n))\\lKS7wo ■Vwo\\L2{n)ds ^ 0, as t +oo, 

since an application of the spectral theorem shows that 

\\e'^^-'^^\\^LHn),LHn)) < sup e~^(*"') ^0, as t ^ +cx). 

AGspec(P) 

This establishes ( 14. 3p . 
Next we shall show that 

\\Pij{t, ■) - -fnVwo ■ Vwo\\L2(^n) ^0, as t ^ +oo. (4.4) 

We have for t > 1, 

P^(t, x) = P f a^(s)e~(*"'^^7KVwo ■ Vw^ds + P [ C'^'^kWwo ■ Vwods. 
Jo Jo 
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The formula 



see [Hj, implies that 
Jo 

and 

||g-(*-i)^'^^yy;p . Vwo||L2(n) 0, as t +00. 

Finally, 



Wo, 







< / a'^{s)\\Pe ''^^\\c{LHn),L^n))hi^'^WQ-VwQ\\L'2(n)ds ^ 0, as t ^ +cx), 



smce 



l|i"e-(*-^)^lU(i2(^),i2(f,)) < sup (Ae-(*"^)") ^ 0, as t ^ +00. 

Aespec(P) 

This proves (14 ■4p and completes the proof of the lemma. □ 
Lemma [4.21 implies that as t — ?■ +00, 

u ■ AV^Plan ^ ■ AViP-\^KVwo ■ Vwo))\an in H'/^dn). 

Thus, as t — 7- +00, we have by a repeated application of the divergence theorem 
together with the Cauchy-Schwarz inequality, 

[ ^^^^,A{a{t)h{x))dS = [ u-AVipdS^ [ u ■ AV{P-\jkVwo ■Vwo))dS 
Jdn ' ' Jdn Jan 

V ■ AV(P"^(7kVwo ■ Vwo))dx = - -fVwo ■ Vw^dx = - / A^{h)hdS. 

Jn Jan 

Here 

h ^ 'jd^wolan, 

is the Dirichlet-to- Neumann map, associated to the problem (14 .11) . Thus, the 
knowledge of the voltage-to-heat flow map ^-y,K.,A determines the Dirichlet-to- 
Neumann map A^. It follows from [38] that the isotropic conductivity 7 G C°°{Q) 
is uniquely determined by S-y,^,^- This completes the proof of Proposition 14.11 
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5. Recovering the heat parameters A and k 
The purpose of this section is to prove the following result. 
Proposition 5.1. = E^ ^^j^^j, then Ki = k,2 o-nd Ai = A2. 

When determining the conductivity 7 in the previous section, we were con- 
cerned with the power densities F which are supported in the region t > 1/2 and 
independent of t near +00. Here we shall instead concentrate the density F in a 
small neighborhood of t = 0. 

Let X e be such that < x < 1, supp (x) C [0, 1] and 

,2 



x\t)dt = 1. 

For e > 0, we define 

Xe{t)=e-'l\{t/e), (5.1) 

so that xi{t) — 7- 5(t), the Dirac measure at t = 0, as £ — )■ 0. Let h G H^{dQ), for 
s large enough. Then the solution u = u'^'^ of the problem 

V-7VM(t,x) = in fi, 

u\dn = Xeit)Kx), 
satisfies u^'^{t,x) = Xe{t)'w^{x), where w^{x) solves 

V ■ jVw'^ix) = in Q, 

y (5-2) 
w Ian = h{x). 

Let h,he H\dVt) and 

pe,h±h ^ ^^^^e,h ^ ^el^ . ^^^e,h ^ _ -^xl{t)^ {w^ ± W^) " ^ {w^ ± w'') . 

Denote by ^^'^^'^^i the solution of ([23D with F = F"''*^'^, k = kj and A = Aj, 
j = 1,2. Set 

a^^\t,x) = a''''''^'^{t,x) = ^{■^''''+^'^{t,x) -ip''''-^'^{t,x)). 

Thus, a^^^ is a solution of the following inhomogeneous initial boundary value 
problem, 

(9, + P>(^-) = i«:,(x)x?(t)7 



= Kj(x)x^(t)7Vw'' ■ Vw'', in M+ x n, 



(5.3) 



"1 \-^JAe 



Since 

P- = -KjV ■ AjV 
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is a self-adjoint positive operator in -^^^ (fi) with the domain Vi^Pj) = {Hq fl 
H'^){Q), the spectrum of Pj is discrete, accumulating at +00, consisting of eigen- 
values of finite multiplicity, < A["''' < Ag"''' < ■ ■ ■ — )• 00. Associated to the 
eigenvalues A^''^ we have the eigenfunctions ip^j^'^ G V{Pj), which form an or- 
thonormal basis in L^^, In what follows we shall assume, as we may, that the 

eigenfunctions (p\^^ are real-valued. Hence, 



k=l 

with convergence in C(M+, P(Pj)). Therefore, 



(3)\\ 
k )\dn 



k=l 



with convergence in H^^'^(d^l), for each fixed t > 1. 
Next we notice that 

1 ~ 1 ~ 

u ■ AjVa^^^lg^^^^ = -J:^^^^^A,{Xs{t){h + h)) - -J:^^^^^A,{Xe{t){h - h)). 



Thus, since S^^^i.Ai = ^j,K2,A2y follows that for all t > 1, 

00 00 

Y^c^lmu . A,V^^^%n = Y.^^lmu-A,V^^;%n. (5.4) 

fc=i fc=i 

Here we would like to let e — 0. In order to do so, it will be convenient to 
obtain an explicit representation of the Fourier coefficients c^^lit). 

Set d\^^ = {'yVw'^ ■ Vw^, y9^-''')^2, where the scalar product is taken in the space 
L'^{Q,dx). It follows from that 

c£(0) = 0. 

Hence, 

4^) (t) = e-^^*4^) r e''^'^xlis)ds, t > 0, 
Jo 

and cj^l{t) is uniformly bounded in k, e. For t > 1 and k = 1,2, . . . , fixed, we get 

c(f ) (t) = e-^i''*4 x\s)e^''^'''ds ^ e-^''d^\ as e ^ 0. (5.5) 
Jo 

Using (15.51) . we may let e — )■ in (15.41) . and conclude that 

00 00 

5^e-^i^M'^(z.-A,V^«)U = ^e-4%l^)(z.-A2V^l^^ t>l. (5.6) 

k=l k=l 
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In what follows we shall have to distinguish the eigenvalues of the operator 
Pj, j = 1,2. In order to do that let us continue to denote by A^"''' the sequence 
of distinct eigenvalues of Pj and let m^"''' denote the multiplicity of A^"'''. Let 
(p\^\, . . . , ^p^''^ be an orthonormal basis of the eigenfunctions, corresponding to 

the eigenvalue A^"'\ By the uniqueness of the Dirichlet series, see 0, and f l5.6p . 
we obtain the following result. 

Proposition 5.2. Assume that S^^^i.yii = ^7,k2,A2- Then for all k = 1,2, ... , we 
have 



X(l) _ x(2) 



and 



(1) (2) 

i). ml 



AiV(^2)|an = E<](^-^2V^S)lan. (5.7) 



i=l i=l 



Here dj:[ = (7VW ■ Ww , (p^^)i2 and h,hE H^(dQ,) are arbitrary functions. 
Let us introduce the following linear continuous operators 



4^' 

Rk\n = J2^F,^^^}),^.{u.A,Vv^^})\an, j = 1,2, k = l,2. 



i=l 



Proposition 13.11 together with (15. 7p implies that i?^^^ = R^^^ on a dense subset of 
L^(fi), and hence, everywhere. On the level of the distribution kernels, we obtain 
that for all A; = 1, 2, ... , 



(1) (2) 



E^2(^)(^ ■ = E^2(^)(^ ■ ^2V^g)(i/),a; G n,ye dn. (5.8) 



1=1 i=l 



We shall next need the following result, see [7]. Since the argument is short, 
for the convenience of the reader, we give it here. 



Lemma 5.3. The functions {u ■ AjViPj^{)\dQ, . . . , (i^ ■ AjV(f j^ (j))\d(^ c-f^ linearly 



independent, j = 1,2. 

Proof. Assume that there are ci, . . . , c (j) G M such that 

U) 



Y,c^{ly■A,Vip)l>)\9n = 0. 



i=l 
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Then = 22i=i ^iVk i satisfies 

PjLp = A^-'V in 

Vlan = 0, z/ ■ AjV(p\dn = 0. 

By the unique continuation principle, we get ip = in Q. Thus, q = for 
1 < z < m^i^^ This proves the lemma. 

□ 



Our next goal is to analyze the consequences of f lS.Sp . and the key step here is 
the following algebraic result which is similar to |71 Lemma 2.3]. 

Lemma 5.4. Let /j : -> M, /j : 9i7 -> M, i = 1, . . . ,m^^\ and gi : Q ^ R, 
gi : on — 7- M, / = 1, . . . , m^'^\ be such that 

= for all X en, y e on. (5.9) 

i=l (=1 

Moreover, assume that the systems {/i, . . . , fmw}, {/i, . . . , fmw}, {gi, ■ ■ ■,gm(2)} 
and {gi, . . . , 'g^(2) } are all linearly independent. Then m^^^ = m^"^^ and there exists 
an m^^^ x m^^^ invertihle matrix T with real entries such that 

F{x) = TG{x), F{y) = {Ty^G{y) for all x en, ye dn. 

Here we use the notation 

f fiix) \ / fiiy) \ 

Fix)= ; ,Fiy)= ^ : 

/ 9iix) \ ^ / giiy) 

G{x)= : ,G{y)= { : 

Proof. As /i is not identically zero in dn, there exists yi e dn such that fiiyi) 7^ 
0. Assuming that 

det { ^^y'} i^yA = for all y e dn, 

we get that /i,/2 are linearly dependent which contradicts the assumptions of 
the proposition. Thus, there exists y2 e dn such that 
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Continuing in the same way, we find points j/i, j/2, . . . , UmW G dQ such that the 
matrix 

n~- •^^'^^2) /2(y2) ••• fmwiy2) 
~ : : : : 

is invertible. It follows from (15. 9p that QjF{x) = QgG{x) for any x G fi, where 

giiVi) ••• 5'm(2)(yi) 



K9i{y 



Thus, F{x) = TG{x) for any x E n, where T = Q~^Qg, and therefore, m^^^ < 

m^'^\ Similarly, using the fact that {gi, ■ ■ ■ ,gm(^)} is linearly independent, we 
have G{x) = TiF{x) for any x E ^l, and m*^^-' < m^^\ Hence, m^^^ = m^'^\ 
Furthermore, 

F{x) = TTiF(x) for all x E Q. (5.10) 

Since the system {fi, ■ ■ ■ , fmw} is linearly independent, in the same way as 
above, we see that there are points xi, . . . , a;„(i) E Q such that the vectors 
F{xi), . . . , F{x^(i)) form a basis in C™'^'. Thus, (I5.10p implies that TTi = / 
and therefore, T is invertible. Similarly, one can see that F{y) = TG{y) for all 
y E dfl with an invertible matrix T. 

It follows from (15. 9 p that F{x) ■ F{y) = G{x) • G{y) and therefore, 

(f *T - I)G{x) ■ G{y) = for al\xEn,yE dn. 

Since there exist points Xi, . . . , E Q and yi, . . . , y^w E dVL such that the 
vectors G'(xi), . . . G(x„(i)) (respectively, G{yi), . . . G{y^(i))) form a basis in C^^^\ 
we get T*T = J. This proves the claim. □ 

If follows from Lemma 15.41 together with (15. 8 p and Lemma 15.31 that m^^-* = 



m 



(2) 



=: rrik, and there exists an rrik x rrik invertible matrix T such that 



X E Q, 





= T 


( ^'Si 


x)\ 






\r k,irn; 


{x)J 



(5.11) 



and 



/z.-Ai(y)V^«(l/)\ 
^VZ/-Ai(y)Vv.l^L,(2/)/ 



/ u-A,{y)V^'-^hy)\ 



{T 



\yA,{y)V^^^(y)J 



y E dn. (5.12) 
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Using (15. lip and fl5.12p . we have 



/z/-Ai(i/)V^g(y)\ 



\u.A,iy)V^^;i(y)/ 



y e on. (5.13) 



The next step is to show that the matrix T is in fact orthogonal. This will 
follow once we establish the following result. 

Proposition 5.5. = S^,k2,A2 then Ai\qq = A2\dn- 

Proof. Consider the following eUiptic boundary value problem, 

PjU = Kj'yVw'' ■ Vw'' in Q, 

u\dn = 0, 



(5.14) 



where (respectively, w^) is the solution to the problem (15. 2p with the boundary 
source h G H'^{dQ) (respectively, h G H^{dQ)), s > 1/2 large enough. We shall 
now return to the original notation, where each eigenvalue A^"'^ of the operator 
Pj is repeated according its multiplicity. Since ^ spec(Pj), the problem (I5.14p 
has the unique solution 

k=l 

with convergence in H^{Q). Thus, 

°° Ai) 

V ■ A,Vu^%^ = Y,^)^ ■ A^^^'i'^W 

k=l ^k 

and therefore, it follows from Proposition 15.21 that if S^^^i.^i = ^7,^2,^25 then 



Define the continuous map 



(5.15) 



do.-, 



U 



is a solution to the problem 



-V • (AjVm(^'^)) = f in 



u 



(FJ) 



an 



It follows from (I5.15P together with Proposition 13. II that "^^^^ = '^^'^^ on a dense 
subset of L'^{Q), and thus, everywhere. Hence, Proposition lA.ll in Appendix A 
implies that Ai\qq = Ailan- 

□ 

Now going back to equation (15.130 . using Lemma [5.31 we obtain that T is an 
orthogonal matrix. Proposition l5.2l together with (15.110 gives the following result. 
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Proposition 5.6. Assume that S^ ^i.^ii = ^^ ^2.^2- Let (/j^^'' be an orthonormal 
basis in L^(fi, dx) of the Dirichlet eigenfunctions of the operator Pi. Then 
the Dirichlet eigenvalues X^^'' (respectively, X^^^ ) of the operator Pi (respectively, 
P2), counted with multiplicities, satisfy X^^'' = Xf^ for all k and there exists 
an orthonormal basis in L'^{^1, K2^dx) of the Dirichlet eigenfunctions (p^^'^ of the 
operator P2 such that ip^g^^ = ipf^ for all k. 

We shall next show that Proposition 15.61 yields that ki = ^2- Indeed, let us 
write 

00 00 00 

k=l k=l k=l 

where the Fourier coefficients Ck are given by 



Cfc 



Kiif'j^^i ^dx = / (p^^^dx = / H2^t^ 1^2 ^dx = dk- 
Jq Jn 



Thus, Hi = K2. It follows that PiU = P2U, for any u G C^{Q), and we get 
Ai = A2. The proof of Theorem 12.11 is complete. 

Theorem 12.21 can be proven by exactly the same arguments presented in this 
section applied to the problem (12.31) with the right hand sides of the form 

F{t,x) = xlmi^), HeL\n), 
where Xe is given by fl5.ll) . 

Appendix A. Boundary reconstruction 

Let Q C M", n > 2, be a bounded domain with boundary, and A{x) = 
{ajk{x)), 1 < j, k < n, he a. real symmetric n x n matrix with aj^fc(x) G C°^(r2). 
Assume that there exists Cq > such that 

n 

('AxMk > Coiei' for all (x, e H x M". 

j,k=i 

Consider the following elliptic boundary value problem, 

-V ■ (AVu) = F in 
u\dn = 0. 



(A.l) 



For any F G the problem f lA.ip has a unique solution u = E H^iVt) fl 

Hq{VI) and one can define the map, 

^ ■ L^in) ^ L^{dn), ^{F) = v AVu^\an, (A.2) 

where v is the unit outer normal to the boundary dVt. We note that the map 
^ is sometimes used to model boundary measurements for optical tomography 
with diffusion approximation, [T71 [HI |35] . 
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We have the following proposition, which is closely related to the earlier bound- 
ary reconstruction results of a Riemannian metric from the Dirichlet-to-Neumann 
map [SI [29]. 

Proposition A.l. The knowledge of the map \1/ given by flA.2p determines the 
values of A on the boundary dQ. 

Proof. We shall recover the values of A on the boundary by analyzing the dis- 
tribution kernel of the map obtained by constructing a right parametrix for 
the boundary value problem flA.ip . Let us denote A = —V ■ (^V). Since A{x) 
is a positive definite matrix, smoothly depending on x, we can view as a Rie- 
mannian manifold with boundary, equipped with the metric G = A~^, G = {gij), 
^ ^ i, j ^ n. To construct a parametrix for ( lA.ll) . we shall work locally near a 
boundary point. Let xq G dQ and introduce the boundary normal coordinates 
y = {y'lUn) ^ U, y' = {yi, . . . ,yn-i), centered at Xq. Here U stands for some 
open neighborhood of in M". In terms of the boundary normal coordinates, 
locally near xq, the boundary dQ is defined by ?/„ = 0, and ?/„ > if and only 
if x e r2. In what follows, we shall write again {x',Xn) for the boundary normal 
coordinates. 

In the boundary normal coordinates, the metric G has the form 

n-1 

G = ^ gap{x)dxadxp + {dxnf, 

a,p=l 

see [29], and the principal symbol of the operator A is given by 

n-1 

a,/3=l 

Therefore, the equation ao(x,^',^„) = 0, ^' = (^i, . . . ,^n-i), has the solutions, 

n-1 

E r^i^)^c.ip. (A.3) 

We can view ^ as a linear continuous map in the space V{U) 
normal coordinates, the problem (lA.ll) has the following form, 

Au = F in M" = {a; G M" : a;„ > 0}, 

m|x„=o = 0. 

Let 

ro(x,0=ao(x,0''(l-x(0), a: G [/, ^ e M", 

where xiO ^ C^(]R"), xiO — 1^1 — 1 X = 1 near 0. The operator 
Op(ro) is a rough parametrix for the operator A, which will be sufficient for our 



In the boundary 
(A.4) 
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purposes. Here we are using the classical quantization of a symbol a G S^{U xM"), 
which is given by 

Op(a)w(x) = -i- / / e'^^-y>^a{x,i)u{y)dydi. 



(27r)" 

As usual, we say that a & S^{U x M") if locally uniformly in a; G f/, we have 

Let r+ be the operation of restriction from to and let e+ be the operation 
of extension by zero from to M". 

We shall construct a right parametrix for the boundary value problem ( ]A.4I) in 
the following form 

i?(F) = r+Op(ro)(e+F) + i?6(^/.). 

Here 

= -ror+Op(ro)(e+F), Tq : m t-> ?^|x'„=o, 
and R}, will be constructed as a right parametrix for the boundary value problem 

Au = {] in W], 

In what follows we shall suppress the operator r+ from the notation, as this will 
cause no confusion. 

When constructing the operator we shall follow the standard approach in 
the theory of elliptic boundary value problems, see [S]. To this end, let x ^ 
C^(M"-i) be such that x = 1 for |f | < 1. Notice that 

(1 - X(r))x(0 = 0. 

Let a = a{x,^') be a simple closed smooth curve in the upper half-plane 
Im^„ > 0, which encircles the root A+(x, in the positive sense. In what 
follows we may and will choose cr so that it is independent of x G f/, depending 
on ^' only, i.e. cr = a{^'). When (p G C^{U (1 R"~^), we define the operator 

(A.6) 

where 

m = [ e-^y'<'ip{y')dy' 

is the Fourier transform of ip. By a contour deformation argument in the complex 
^„-plane, we have 

liip = -Op((l - x)%^){v ® 5x'„=o), Xn > 0. 
t 

We get therefore, 

AUif = Op(6)(v9 ® 4„=o), b G S-\U X R"), xn > 0, 
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since the operator A is local. 

We shall take Rb{ip) = Hyj, for some function ip, defined locally near G M""^, 
to be found from the boundary condition, i.e. 

ToU^ = ij. (A.7) 

To this end, we shall now prove that roll is an elliptic pseudodifferential operator 
on the boundary and compute its principal symbol. By the residue calculus, 
using (]A.6p . we get 

r„n = 0p(4). M^,e) = (1 - M^)"^" 

where 

n-l 

We introduce next a rough parametrix of roll, given by Op{do), where do G 
S\{Un M"-^) X R'^-i) is such that 

do = 2i\^'\A, for |^'| large. 

To satisfy ( ]A.7p . we choose 

<f = Op{do)^ = -Op(rfo)(roOp(ro)(e+F)). 

This choice of ip completes the construction of a rough parametrix for the bound- 
ary value problem ( 1A.5I) . given by 

Rb{iP) = -nOp(^o)(roOp(ro)(e+F)). 

Hence, the parametrix for the problem ( 1A.4I) has the form 

R{F) = Op(ro)(e+F) - nOp(rfo)(roOp(ro)(e+F)). 

In the boundary normal coordinates, the operator \E' is given by 

and therefore, to obtain the claim of the proposition it suffices to analyze the 
distribution kernel K{x',y) of the operator Todx„R given by 

[TodxMF)){x')= f K{x',y)F{y)dy, x'eW-\ yER^. 

Let us first consider the Schwartz kernel of the operator 

ro5.„Op(ro)(e+F) = -i-ro9,, / e^^Vo(x, Oe^^lO^^, 
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which is given by 

Recall that here |/„ > 0. Restricting the attention to the region |^'| > 1, by a 
contour deformation argument to the lower half plane, we find that 

/ e-^"«"ie„ro(x', 0, O^Cn = 2me-'y-^- , f ^ , 

and therefore, 



p-ynW\A / 1 \ 

t{x\ Vn, a = + e-^"l« '"^ j ' lel > 1- 



Next, the operator rodx„U. is a pseudodifferential operator on ^, given by 

The principal symbol of the operator ro(?a;„n is therefore 

^ ^ rc^^n = , ^, , , = J, I'^l large enough. 



27ri ao(a;',0,0 a^„ao(a;', 0, C', A+) 2' 

The operator Op((io) is also a pseudodifferential operator on ]R"~^ and its princi- 
pal symbol is given by 2A+ G S'^(([/ flM"^"'^) xR"^~^), |^'| large enough. Hence, the 

principal symbol of the operator roc?a;„nOp((io) is iA+ G S^{{U fl M""^) x M""^), 
and therefore, its kernel is given by 

e') = iA+(x', 0, + Pq{x', e'), IC'I large enough, 

where po e n R"-^) x R"-^). 

Finally, the kernel of the operator ToOp(ro) is given by 

K?\z\y) = j:^ [ e'^^'-y'>^'e-'y-^-ro{z',0,v)d7j 

\ Q-y-nW\A 

d2{z', Vn, v') = o—TT^ > IVI large enough. 
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Here as usual we use the residue calculus, where only the pole in the lower half 
plane contributes. 

Hence, the kernel of the composition ro(9a;„HOp((io)ToOp(ro)(e+F) is given by 
K^'\x\y) = I Kl'\x',z')Ki'\z',y)dz' 
1 



(27r)2("-i) 
1 




eii^'-'')<'e'^''-y'>'^'di{x',C')d2{z',yn,v')d^'dv'dz' 



where 

c(a;',y„,r^')= j j e^^-'-'''y^^'-^'^d,{x\0d2{z\yn.v')di'dz' . 

Here ?/„ > occurs as a parameter. Now di{x' G S^{{U fl M"^^) x M"-"^) for 
large |^'| and d2{z' ,yn,ri') satisfies 

\d:,d^^,dld2iz\yr,,r^')\ < Q,,,e-^"<'''>(r^r 



where {t]') = ^/T+~\Tf\^ is large enough. Hence, c{x' ,yn,rj'), depending on the 
parameter is the symbol of the composition of two pseudodifferential opera- 
tors in the tangential directions. By the standard results on pseudodifferential 
operators, see [H], it has the following asymptotic expansion, 

c{x',yn,v') ~ Y] ^d^'di{x',ri')D^,d2ix',yn,v'), 

\a\>0 

with the leading term di{x' ,ri')d2{x' ,yn,f]')- The knowledge of the operator \E' 
implies the knowledge of the kernel K^^\x', y) — K^'^\x', y) for any x' G t/nM"^^ 
and y eU r\ M" . This implies the knowledge of 

_l_ g-j/n|5 \aq( j _ c{x' , yn, C,') , for any |^'| large enough. 

Vis U/ 

The leading term of the latter expression is given by 

L _ p-J/n|5'U 

2 2 

Varying we recover A{x', 0). The proof is complete. 

□ 
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